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ABSTRACT 

We report on a signature of chromospheric downflows in two emerging-flux regions detected by 
time-distance helioseismology analysis. We use both chromospheric intensity oscillation data in 
the Ca II H line and photospheric Dopplergrams in the Fe I 557. 6nm line obtained by Hinode /SOT 
for our analyses. By cross-correlating the Ca II oscillation signals, we have detected a travel-time 
anomaly in the plage regions; outward travel times are shorter than inward travel times by 0.5 — 1 
minute. However, such an anomaly is absent in the Fe I data. These results can be interpreted 
as evidence of downflows in the lower chromosphere. The downflow speed is estimated to be 
below 10 km s _1 . This result demonstrates a new possibility of studying chromospheric flows by 
time-distance analysis. 

Subject headings: Sun: chromosphere — Sun: faculae, plages — Sun: helioseismology 



1. INTRODUCTION 

Local helioseismological techniques have been exploited to reveal subsurface structures and dynamics 



of the Sun, for instance, mass flows and wave-speed variations below sunspots (e.g., Zhao et al. 20011, or 



subsurface structure of supergranular flows (e.g., [Sekii et al.j|2007 l. In this Letter, we examined whether we 



can detect flows in emerging magnetic flux regions. 

Magnetic flux tubes are generated in the convection zone. They emerge into the photosphere due to 
magnetic buoyancy, forming sunspots and/or restructuring the coronal magnetic configuration, eventually 
leading to activity phenomena, such as flares. Emerging flux regions are the place where such flux tubes 
are first observed directly; their configurations and dynamics are manifestations of interaction between 
plasma and magnetic flux in the subphotosphere. Thus, studying various aspects of flux emergence is vitally 
important in understanding solar and stellar dynamo mechanisms. 

Numerous studies about emerging flux regions have been done both observationally and theoretically, 
and revealed complicated structures and flow patterns (e.g., Strous et al.|1996 Kosovichev et al.|2000 Kozu 
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et al.||2006| |Pevtsov fe Lamb||2006l |Magara|[2008l |Cheung et aT1|2008|). It is generally thought that along 



the flux tube are plasma downflows. The evacuated plasma decelerates as it descends due to the increasing 
density; in fact, the observed speed is 30 — 50 km s _1 in the chromosphere and 1 — 2 km s _1 in the photosphere 



(e.g., Tajima fe Shibata||2002 and references therein). To understand how emerging flux regions evolve and 
how they affect activity phenomena, it is essential to study them in a range of layers from the subphotosphere 
to the upper atmosphere. 

Local helioseismology is the only way to investigate the subphotospheric structure and dynamics of the 
Sun. Moreover, it allows us to study the dynamics above the visible surface. In this Letter, we report 
on a signature of chromospheric flows detected by a time-distance helioseismology analysis. Examining 
two emerging flux regions observed by Hinode, we detect a signal in the travel-time differences, which may 
be caused by plasma downflows in the chromosphere. Another implication of the present work is that 
chromospheric dynamics must be taken into account when chromospheric lines are used in helioseismology 
analysis. We briefly introduce Hinode observations in Section [2] and describe the analysis methods and 
results in Section [3] Discussion of our interpretation and conclusion are given in Section [4] 



OBSERVATIONS 



The Solar Optical Telescope (SOT; Tsuneta et al. 2008) onboard the Hinode satellite (Kosugi et al. 



2007) made two observing runs for helioseismology around the solar disc centre: on 2007 November 23 
(11:17-23:55 UT) and on 2007 December 04 (0:14-23:58 UT). Hereafter, we refer to the first 12-hour run 
as data set 1, and the second 24-hour run as data set 2. In the field of view of data set 1, an emerging 
flux region appeared before the observation started, showing a plage and small sunspots. This region was 
registered as NOAA Active Region 10975. It disappeared within two days after the Hinode observations. In 
the field of view of data set 2, a decaying active region (NOAA Active Region 10976) was observed. Sample 
images of these regions are shown in Figure [I] 

During these observations, both SOT instruments, the Narrowband Filter Imager (NFI) and the Broad- 
band Filter Imager (BFI), were used: the BFI took Ca II H images, while the NFI took nonmagnetic Fe I 
557.6 nm filtergrams. The time cadence was 60 seconds. The Fe I Dopplergrams are calculated from the 
blue-wing and the red-wing intensities of the Fe I 557.6 nm line; that is, the (blue - red) signal is divided by 
the (blue + red) signal. During the each observing period Hinode tracked the regions using the correlation 
tracker (CT; |Shimizu et"aT1|2008| ) to stabilize the images. The field of view is 218 x 109 arcsec for the Ca 



II H data and 328 x 164 arcsec for the Fe I data, while the pixel sizes are 0.108 arcsec (for Ca n H) and 
0.16 arcsec (for Fe i). In this analysis, we use only the field of view of the Ca II H data. The Fe I data are 
interpolated to the Ca II H data grid. We also apply 2x2 pixel summing, increasing the pixel scale to about 
0.2 arcsec. For the intensity data sets, we use running difference of images to remove possible spatial trend. 



TIME DISTANCE ANALYSIS AND ITS RESULTS 



We analyze the 12-hour and 24-hour data sets (data sets 1 and 2) by a time-distance helioseismology 
technique ( puvall et al. 19931. First, we apply three filters in the frequency-wavenumber space: (1) an 
f-mode filter to remove the f-mode signal, since we use p-mode waves only, (2) low-wavenumber filters and 
a Gaussian frequency filter with a 1-mHz width peaking at 3.3 mHz to remove the known artifacts (see 
Tsuneta et al.|[2008 1 and select the main p-mode signal, (3) a phase-speed filter with a width of 4.4 km s _1 
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centered at a speed of 22 km s 1 to select waves travelling to a distance of 14.86 Mm, which is chosen for 
this study. 

Second, by cross-correlating the filtered signals, we measure the travel-time difference between the out- 



ward and inward wave components following the standard time-distance helioseismology procedure (Koso- 



vichev fe Duvall|1997 1. For this, we average the signals over an annulus of the radius of 14.86 Mm around a 
target point, and cross-correlate the signals at the target point and that of the annulus. We measure outward 
(from the target to annulus) and inward (from annulus to the target) travel times at each point in the field 
of view by fitting a Gabor- wavelet function to the cross-correlation function. The fitting function, C(A,r), 
for spatial displacement A, and time lag r has the form 



C(A, r) = Acos(o;o(t — r p )) exp 



5u> 



(1) 



where A is the amplitude, loq is the central frequency, Su is the frequency bandwidth, and r g and t p are the 
group and phase travel times. These five parameters are determined by applying a nonlinear least-squares 
fitting method (Levenberg-Marquardt method; see, e.g., Press et al.||1992 ). Errors of the fitted parameters 
are derived from the estimated covariance matrix. It should be mentioned that we used the phase travel 
time in the following discussions, because the phase travel time is determined with higher precision than the 
group travel time. 

Figure |2]shows outward-inward travel-time difference maps. In these maps, misfitted points are assigned 
null travel-time differences; the number of misfittings is insignificant. In the quiet regions, both the Ca II 
H intensity data and Fe I Doppler data show the same pattern of flow divergence due to supergranulation. 
This is similar to the results from a previous Hinode observation reported by Sekii et al. (20071. In the 



plage region, we find a strong travel-time anomaly in the Ca II H data: the outward travel time is shorter 
than the inward travel time by about 0.5 - 1 minute on average. This anomaly is not present in the Fe I 
557.6 nm Doppler data. We also analyze the Fe I 557.6 nm intensity data from the same observing runs. 
These data are noisier, but they show the travel-time difference maps consistent with these obtained from 
the Fe I 557.6 nm Dopplergrams. Particularly, they do not show the travel-time anomaly, cither. 

The anomaly in the outward-inward travel-time difference in the chromospheric data is seen also in the 
averaged cross-correlation function. We average the cross-correlation functions in the quiet Sun (QS) region 
and in the plage for both data sets 1 and 2. We define the plage as a region where the Ca II H intensity is 
25% larger than the average intensity of QS. By fitting the Gabor wavelet to the averaged cross-correlation 
functions, we obtain the travel times listed in Table [T| We find that in the plage the outward travel time is 
about 0.67 min (for data set 1) or 0.42 min (for data set 2) shorter than the inward travel time only in the 
chromospheric Ca II H data, while in the QS the difference between the outward and inward travel times is 
insignificant in both the chromospheric and the photospheric data. 

We also have checked temporal variation of the travel-time anomaly for data set 2. The 24-hr run is 
divided into three 8-hr runs, and the travel times arc measured for each 8-hr run. As a result, we obtain 
the following sequences of the outward-inward travel-time difference in the chromospheric Ca II H data: 
(-5.35 ± 0.04) x 10" 1 min, (-4.16 ± 0.04) x 10" 1 min, and (-4.21 ± 0.06) x 10" 1 min in the plage region, 
and (-1.11 ±0.03) x lO" 1 min, (-0.14 ±0.07) x lO" 1 min, and (+1.37 ±0.03) x lO" 1 min in the QS. In the 
plage, for all the 8-hr runs the anomaly similar to the one detected in the entire run is detected again, and 
particularly it is rather stronger in the first 8-hr run. This variation is consistent with a decaying behaviour 
of the plage region if we interpret that the travel-time anomaly results from a chromospheric downflow 
(see Section [4] for details), and that the downflow weakened during the observation. We also find, in the 
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developing plage in data set 1, the anomaly being stronger in the latter half of the observation. In the QS, 
the weak travel-time anomaly simply corresponds to the supergranular pattern. The observed variation with 
time implies that the lifetime of supergranules is slightly less than 1 day. This is consistent with previous 
studies (e.g.,|Stix||2002|). 



4. DISCUSSION AND CONCLUSION 

In our analysis, we have used only acoustic (p-mode) waves with the central frequency of 3.3 mHz. 
Although stationary harmonic p-mode waves at this frequency are evanescent in the chromosphere, in reality 



such perturbations excited by impulsive subphotosphere sources do propagate upward (e.g., Lamb 1975 
Goode et al"1|1992| ). This certainly requires further careful modelling, but here for the initial interpretation 
we assume that the waves travel roughly with the chromospheric sound speed. Thus, a chromospheric signal 
is observed after the corresponding photosphcric signal is observed and the time delay is the acoustic travel 
time between the two layers, which is estimated by L/c ~ 25 second, where L ~ 250 km is the formation 
height difference between the Fe I 557. 6nm line and the Ca II H line ( Carlsson et al.|2007 l and c ~ lOkms -1 



is the sound speed in the chromosphere. Indeed, the phase shift between the photosphcric G-band intensity 
signal and chromospheric Ca n intensity signal around 4 mHz on the p-mode ridge, reported by |Mitra-Kraev 



et al. (20081, implies a ~ 30-second propagation time between the photosphere and the chromosphere. Such 
a direct measurement in our case is not straightforward because our photospheric signal is Doppler velocity. 
Understanding phase shifts between Doppler signal and intensity signal would be part of a future work. 
Please note, however, the travel-time analysis we use is not affected by such a difficulty as we cross-correlate 
the same observables. 

Figure [3] shows how the disturbance caused by a source at point C propagates. If the wave arrival times 
at points A and C in the photosphere are t a and tc , the travel time from C to A is simply tc to A = t a — tc ■ 
The 'travel time' from C' to A', defined through the cross-covariance fitting procedure, tctoA'; is identical 
to tctoA if the time delays at A, dt\, and at C, die, are the same. Our travel-time measurements in the 
QS are consistent with this picture; the photospheric and the chromospheric travel times arc essentially the 
same (see Table [j}. 

If there is a downflow above C, however, and there is no downflow above A, then the time delays at two 
points are different: dtc > dtA- By setting point C as the center of an annulus, or as the target point, and 
point A as a point on the annulus around the target point, we can describe the outward travel time at the 
photospheric level, t ut.phj and the travel time at the chromospheric level, i ou t, c h> as 

tout.ph = ictoA = ^A — tc (2) 
tout.ch = tC'toA' = tCtoA + {dtx — dtc) < tctoA j (3) 

and also the corresponding inward travel times ti n , P h and ti n ,ch as 

tin.ph = ^AtoC = tctoA (4) 
iin,ch = *A' toC = ^AtoC — (^A — dtc) > ^AtoC j (5) 

where we assume for simplicity that there is no perturbation in the subphotosphcric layer, or, cquivalcntly, 
removed the subphotospheric components from the discussion. Then, the outward-inward travel-time dif- 
ference in the chromosphere, At c h, is 



At c h = t out .ch — ijn.ch — 2(dtA — dtc) < 0. 



(6) 
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The time delay difference between the two points dt A — dtc can be estimated by 

L V 

dt A -dt c = , (7) 

c c 

where V is the downflow speed. Using Equations (J6j) and ([7]), we can roughly estimate the downflow speed: 
at V ~ 8 km s _1 for the plage in data set 1 (Ai c h ~ —0.67 minute) and at V ~ 5 km s -1 for the plage in 
data set 2 (Ai c h ~ —0.42 minute). 

Note that Equatios ^ and ^ imply that the chromospheric mean travel time (tout.cb + £in,ch)/2 should 
essentially be identical to the photospheric mean travel time (t ou t, p h + ^in,ph)/2- Table [I] shows that this is 
indeed the case. 

Actually, we do detect the photospheric downflows in the plage region in the photospheric Fe I Doppler- 
gram. Presumably these flows are much weaker than the flows in the chromosphere. Although we still need 
a more fine-tuned calibration of the NFI data to derive the better estimates of the speed, a redshift signal 
is clearly seen in the Dopplergram (see Figure [T| . Since the signal is somewhat comparable though weaker 
than the Doppler shift due to the solar rotation, the downflow speed is roughly estimated at a few hundred 
m s . This is consistent with the deceleration picture in emerging flux regions. 

In our analysis, we use 14.86-Mm annulus to measure the inward and outward travel times. We chose 
this size because we know that this size of annulus is suitable for detecting spatial patterns of the size of 



the supergranular flow, as is mentioned in the previous works by Sekii et al. (20071. We also check the cases 
for different sizes of annulus for data set 1: from 4 Mm to around 20 Mm. Phase-speed filters are changed 
depending on the annulus size. For these calculations, so that we can measure travel times for smaller annuli, 
we use a Gaussian frequency filter 2.4-times broader than that mentioned in Section 3, although it renders 
wavepackets narrower, and makes the phase travel-time measurement less precise. For the case of annuli 
larger than the original one, the low-wavenumber filter, required for removing the artifacts in Fe I data, 
cannot be used. Thus, for these annuli the travel-time can be measured only for the Ca II data. Results 
are shown in Figure [4] In the QS, the travel-time difference remains small regardless of the size, while the 
difference in plage increases as the annulus size increases. This can be explained by increasingly large part of 
an annulus located outside plage where there is no downflow, which contributes to the travel-time anomaly. 
In fact, when the size attains the spatial extent of the plage (~ 15 — 20 Mm), the anomaly seems to reach a 
plateau. 

We note that a sound speed anomaly in the chromosphere can also lead to an outward-inward travel- 
time difference through affecting the chromospheric time delay. Such an anomaly may arise from magnetic 
field. However, generally in plage magnetic field is not strong enough to explain the travel-time difference 
we have found. Indeed, in MDI ( Scherrer et al. |[T"995 1 magnetogram the plage magnetic field in data set 1 is 



a few hundred Gauss, which translates to only a few hundred m s _1 (magneto-)sound-speed anomaly, and in 
data set 2 it is even weaker. A temperature increase in the plage does not explain the travel-time difference 
cither, since it will decrease the time delay rather than increase it. 

Let us also consider the effect of horizontal motions in the plage. The plage regions changed their shape 
and the bright patches and/or the small sunspots in the field of view were moving during the observation 
periods. The speed of these motions was ~ 1 km s _1 at most, and this is slower than the sound speed by 
an order of magnitude. Our estimates show that these motions are too slow to affect the estimates of the 
chromospheric downflows in any way. 



In this study, we have found a signature of chromospheric downflows in the differences between the 
acoustic travel times measured simultaneously in the photosphere and chromosphere. The outward-inward 
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travel-time anomaly is found only in the chromospheric measurements of a plage region, and, therefore, this 
can be interpreted as chromospheric downflows in the plage. This kind of signature reminds us that chromo- 
spheric helioseismology data may include not only information about subsurface structures but also about 
chromospheric structures and dynamics. This needs to be taken into account in analysis and interpretation 
of helioseismic observations in the chromosphere. Such effects have not been sufficiently studied. 

On the other hand, this result opens new possibilities of multiwavelength time-distance helioseismology 
studies of chromospheric flows and structures. As recent high-resolution observations have revealed, the 



chromosphere is full of activities and flows (e.g., Katsukawa et al. 20071. Multiwavelength helioseismic 



observations may provide a new method for detecting chromospheric flows, which may be of great importance 
for revealing dynamics of the chromosphere. 
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Fig. 1. — Sample images from two data sets. Top row: plage region observed on 2007 November 23 (data 
set 1), Ca II H line intensity image (a) and Fe I 557.6 nm Dopplergram (b). Bottom row: eecaying region 
observed on 2007 December 04 (data set 2), Ca II H intensity image (c) and Fe I Dopplergram (d). In the 
Dopplergrams, white means a redshifted signal, while black means a blueshifted signal. 
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Fig. 2. — Outward-inward travel-time difference maps. Top row: data set 1, (a) Ca II H line intensity and 
(b) Fe I 557.6 nm (Dopplergram). Bottom row: data set 2, (c) Ca II H line intensity and (d) Fe I 557.6 nm 
(Dopplergram). The gray scale corresponds to the range from —1 minute to +1 minute. The field of view is 
the same as in Figure [T| 
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chromosphere 




Fig. 3. — Schematic illustration showing wave propagation (ray path) from the subphotosphere to the 
chromosphere. If the wave arrival times at A and C in the photosphere are t\ and tc, the disturbance 
reaches at A' and C in the chromosphere after some time delay, i.e., at t\ + dt\ and tc + dtc (see the text 
in Section [l]). 
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Fig. 4. — Outward-inward travel-time difference for data set 1 plotted against annulus size. The crosses 
and asterisks indicate the difference measured in plage and QS, respectively. The gray solid lines indicate 
the differences for Fe I Dopplergram. The dotted and dashed lines are for the Ca n intensity data; the former 
and the latter are the differences measured with and without low-wavenumber filters (see the text for the 
details). Original values arc also shown at 14.86 Mm, which are measured using a frequency filter different 
from others, although the values with different filters are essentially the same. Error bars are smaller than 
the symbol size. 



Table 1. Travel Times (in minutes) for Various Data Sets, Regions and Spectral Lines 





Data Set 1 


Data Set 2 




Ca II H 


Fe I 


Ca II H 


Fe I 


Plage 

In 
Out 
Out-In 


22.043 ± 2.8E-03 
21.372 ± 3.0E-03 
-6.72E-01 ± 4.1E-03 


21.725 ± 1.8E-03 
21.646 ± 2.0E-03 
-7.82E-02 ± 2.7E-03 


22.104 ± 3.0E-03 
21.685 ± 3.4E-03 
-4.19E-01 ± 4.6E-03 


21.788 ± 1.8E-03 
21.892 ± 2.0E-03 
1.03E-01 ± 2.7E-03 


QS 

In 
Out 
Out-In 


21.964 ± 2.1E-03 
22.027 ± 3.3E-03 
6.34E-02 ± 3.9E-03 


21.806 ± 1.5E-03 
21.939 ± 2.1E-03 
1.33E-01 ± 2.6E-03 


22.232 ± 2.8E-03 
22.249 ± 3.2E-03 
1.66E-02 ± 4.3E-03 


22.019 ± 2.2E-03 
21.993 ± 2.1E-03 
-2.55E-02 ± 3.1E-03 



